On equilibration and coarsening in the quantum O(N) model at infinite N 
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The quantum O(N) model in the infinite N limit is a paradigm for symmetry-breaking. Qualita- 
tively, its phase diagram is an excellent guide to the equilibrium physics for more realistic values of 
N in varying spatial dimensions (d > 1). Here we investigate the physics of this model out of equi- 
librium, specifically its response to global quenches starting in the disordered phase. If the model 
were to exhibit equilibration, the late time state could be inferred from the finite temperature phase 
diagram. In the infinite N limit, we show that not only does the model not lead to equilibration on 
account of an infinite number of conserved quantities, it also does not relax to a generalized Gibbs 
ensemble consistent with these conserved quantities. Nevertheless, we still find that the late time 
states following quenches bear strong signatures of the equilibrium phase diagram. Notably, we find 
that the model exhibits coarsening to a non-equilibrium critical state only in dimensions d > 2, that 
is, if the equilibrium phase diagram contains an ordered phase at non-zero temperatures. 



I. INTRODUCTION 

The study of the out of equilibrium dynamics of closed 
quantum systems has intensified greatly in recent years 
inspired in considerable measure by advances in the ex- 
perimental study of cold atomic system*^ (see Ref. |BHH1 
for recent reviews) . One important theme in this work is 
the presence and nature of equilibration starting from a 
non-equilibrium state and its relationship to integrabil- 
ity. An interesting stream of work has postulated and 
examined the notion of a generalized Gibbs ensemble 
(GGE) in which an integrable system relaxes to a max- 
imum en tropy state consistent with all of its constants 
of motiorPES) A second important theme is the inter- 
play between non-equilibrium dynamics and phase struc- 
ture w hich g oes under the moniker of the Kibble-Zurck 
probleni 11 ! 12 !, wherein the system is driven out of equilib- 
rium by non-adiabatic changes of parameters, typically 
between values in different phases or between a phase 
and a critical point, see e.g. Refs. IT3"r[T31 

A large fraction of this work has concerned itself, for 
natural reasons of tractability, to systems in spatial di- 
mens ion d = 1 where analyticP*"^]] anc j computational 
, n Kg can j^ readily brought to bear. This restric- 
= 1 however does not allow the study of the 



power 
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impact of dimensionality, known to be important for 
equilibrium behavior. In this paper, we study a model 
where one can move between dimensions while retain- 
ing tractability — the quantum O(N) vector model in the 
infinite N limit (see Ref. [2D for a comprehensive intro- 
duction). This model has much to commend it. It yields 
an equilibrium phase diagram in the Gibbs ensemble in 
various dimensions whose topology is correct for N > 3. 
It also yields critical exponents which incorporate cor- 
rections beyond Gaussian critical behavior in d < 3 and 
correctly locates the lower and upper critical dimension 
at d = 1 and d = 3 respectively. 

The above results assume that the Gibbs ensemble is 
reached. In this paper we study the late time states of the 
model starting out of equilibrium and ask to what extent 



the equilibrium phase diagram is a guide to the late time 
behavior. We note that the model has certainly been the 
object of prior study, in the first instance from a cosmol- 
ogy inspired interest in non-equilibrium field theor y 25 " 28 ^ 
and more recently from the condensed matter/statistical 
mechanical viewpoinlPail]] We build on this work but 
find that there are still new things to say on this prob- 
lem, largely due to asking some fresh questions from the 
condensed matter/statistical mechanical perspective. 

Our results organize themselves naturally into the two 
themes we noted at the outset. First, we revisit the ques- 
tion of equilibration in the infinite N vector model. It has 
been observed previously that the model has a large num- 
ber of conserved quantities that will prevent a relaxation 
to a Gibbs ensemble. We show further that the model 
also does not relax to the GGE that builds in these con- 
served quantities. This behavior is in contrast to that of 
the non-interacting, purely Gaussian, vector model which 
does relax to a GGE , thou gh now with a larger number of 
conserved quantifie d 32 * 33 !. Second, we examine spatially 
homogenous quenches starting from ground states in the 
disordered phase to the critical coupling and into the cou- 
plings in the ordered phase. For the latter, we show that 
the late time state exhibits coarsening in the sense of a 
diverging equal time correlation length in d > 2 for suf- 
ficiently gentle quenches. Our numerics further suggest 
that the system coarsens towards a non-equilibrium crit- 
ical state as t —>■ oo. For larger quenches in d > 2 or 
for any quenches in d < 2, we find that no coarsening is 
possible. Instead, the late time state is disordered. For 
quenches to the critical point, we find no coarsening for 
all d but coarsening in the scaling limit for d > 3 con- 
sistent with the lack of scattering at the Gaussian fixed 
point and previous results on the Gaussian theory^. In- 
terestingly, the results on quenches are qualitatively what 
one would predict assuming equilibration following an in- 
jection of energy density and yet they hold for a system 
that does not equilibrate. 

In the following, we document these claims. We be- 
gin with a review of the infinite N vector model and the 
equations that must be solved to determine its dynamics 
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FIG. 1. Topology of the equilibrium phase diagram of 
the O(N) model in the infinite N limit in spatial dimen- 
sions 1 < d < 2. The dashed lines indicate the different 
quenches that we study. 



FIG. 2. Topology of the equilibrium phase diagram of the 
O(N) model in the infinite TV limit in spatial dimensions 
d > 2. The dashed lines indicate the different quenches 
that we study. 



in Section [XTJ In Section |III[ we discuss the issue of er- 
godicity or equilibration. In Section [IV| we describe our 
results on global quenches into the ordered phase. In Sec- 
tion |Vj we discuss global quenches to the critical point 
and show that the stability of the Gaussian fixed point for 
d > 4 when it comes to equilibrium behavior also shows 
up in the behavior of non-equilibrium quenches provided 
one takes a scaling limit discussed in Ref [2J We finish 
with some concluding remarks in Section [VT| and relegate 
some technical material to an appendix. 



II. THE MODEL 

The Hamiltonian of the quantum O(N) model in d 
spatial dimensions is 



\Jd d x (jll| 2 + |V$| 2 + r|$| 2 + A|$|4^ ( (1) 



H =2 



where $ and II are canonically conjugate ./V-component 
fields, 



[$ i {x),Il j (x')]=i6 d (x-x')6 i: 



(2) 



In the limit N — > oo, the equilibrium physics is soluble. 
In the disordered phase, ($) = 0, the ground state of H is 
well-approximated by the ground state of a free field the- 
ory with mass m c ff that is determined self-consistently: 



(3) 

(4) 



where A is the cutoff. The effective mass m e s is propor- 
tional to the single-particle gap and controls the phase 
diagram (Figs. [I] [2]). In the disordered phase, m e s > 0, 




while in the ordered phase, when suitably re-defined, m e g 
is the mass of the Goldstone bosons and equals zero. We 
describe the phase diagram in greater detail in Sec. |IV| 
The solubility in equilibrium can be traced to two related 
sources: 1) The ground state is a (Gaussian) symmetric 
product state over component indices, and 2) the expec- 
tation value: 



(5) 



factorizes to leading order in 1/N. These two features of 
the infinite N model in equilibrium make the dynamical 
problem tractable as well. 

Now consider preparing the system at t = in a prod- 
uct state symmetric in the component indices^ 




JV 



iv-(o)) =ni 



(6) 



i=\ 



The wavefunction of each component, \j3), can be freely 
chosen. In particular, |"0(O)) need not be Gaussian or be 
a ground state of the infinite N model anywhere in the 
disordered phase. In general, it does not satisfy Wick's 
theorem. The evolution of the state \i/j(0)) for t > is 
generated by H. The Heisenberg equations of motion 
are: 



dt 
dt 



IT,; 



V 2 $, - r$j 



iV 



|$M$,; 



(7) 



For the rest of this article, we work in the Heisenberg 
picture. Remarkably, the factorization in equilibrium in 
Eq. (l5| holds out-of-equilibrium as well: 
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(8) 



As in equilibrium, this factorization leads to an effec- 
tive mass in a free field theory that is determined self- 
consistently. However, the effective mass is now time- 
dependent. More formally, consider the Hamiltonian for 
a free field theory with the time-dependent mass m c g(£) 
determined self-consistently at each t: 



H eS (t) 



d d k 

(2n) d 



l%| 2 



"&(*)) I *jfl 



m 2 ee (t) = r + X' 



N 



(9) 
(10) 



Above, [^ l (k),U j (k / )} = i(2n) d 5 d (k- k')8 VJ . We evolve 
the state \ipifi)) with H c g. We also evolve 1^(0)) with H 
in the limit N — >• oo. The formal statement is that all 
correlation functions involving a finite number of com- 
ponents are identical in the two cases at any fixed t. 
Thus, to determine observables and correlation functions, 
we need only solve for the dynamics in a free field the- 
ory with mass m c ff(£). As in equilibrium, the infinite 
N model out-of-cquilibrium goes beyond the free field 
theory through the single self-consistency condition on 
m G ff(t) (Eq. (lOl). Henceforth, we suppress the compo- 



nent index when an expectation value is independent of 
it. Specifically, we replace (\<&(t)\ 2 )/N by the expectation 
(<& 2 (£)) of a single component. When not indicated, all 
operators in a correlation function are assumed to have 
the same component label. 

A convenient way to determine correlation functions 
from H c s(t) is to expand <&z(t) and Ils(t) in a fixed basis 
at t = 0: 



*s(*) 



%(*) = 



dt 



W) t 

=— a' - 

V2 - k 



(11) 



Above, the (a^,a-) and (a_^,a -) are two independent 
sets of fixed ladder operators with the usual commuta- 
tion relations. The fz(t) are complex- valued coefficients 
known as mode functions. The commutation relations 
stipulate that: 



Im[/g(t)/|(t)] = l. 



(12) 



Given an initial state, we may choose any fixed basis to 
decompose the field operators in. This is a coordinate 
choice; in this article, we pick the particular fixed basis 
in which: 



(<W> = o, /g(0) = /2(0). 



(13) 



It is then easy to see that ft(Q) and |/s(0)| fix all the 
coefficients in Eq. (Ill at t = 0. For example, Eq. (12) 



determines the phase of /r(0) etc. Both functions follow 
from the two-point functions in the initial state: 

<%(0)4,(0)) - JkK L k (2n) d S d (k-k>) 



■/.>'V"jj,w/- 2 



(%(0)Hl(0)) = 



{2TT) d S d (k - k') 

He[<$g(0)%(0))] = Re[f R (0)}f % (0)^(2n) d 8 d (k-k'). 

(14) 



if we substitute for the phase of /e(0) from Eq. ( 12 ). A/V 

above is related to the sum of the occupations of the ±fe 
modes in the initial state: 

Af u =l + n u + n_ % , (a\ H ,) = n~ k (2n) d S d (k k'). 



Finally, from H c ^(t), the definition of m 2 s in Eq. (10 1, 
and Eq. ( 11 ) , we obtain the equations of motion for the 



mode functions: 



dt 2 



+ |*r + mk(*)/jK*)=° ( 15 ) 



™ 2 cs(t)=r + 



X 



d d k 

J2^y 



li^-% 



(16) 



Given the correlation functions in Eq. ( 14 ) and the con 



dition in Eq. (12), we can propagate this system of cou- 
pled second-order differential equations forward in time 
to solve for fj:{t). We may then compute any correla- 
tion function of interest by means of Eq. (11) and the 
solutions for the mode functions. 



III. ERGODICITY 

To determine the nature of the late time state, we need 
to identify all the quantities that are conserved in the dy- 
namics (Eq. ( 15|16 >). First, as the evolution is generated 
by a Hamiltonian, the total energy is conserved. This 
would be true for any isolated quantum system. Addi- 
tionally at infinite N, the operators LS defined below also 



commute with H in the class of states in Eq. (|6| . 



L l 



$tnt - $,-n r 

k k k k 



Thus, 



dt 



0. 



(17) 



(18) 



We use k > to denote half the momenta; each k in the 
set labels the pair (k, —k). The conserved quantity (L|) 
has two interpretations. One is the angular momentum of 
the 2d harmonic oscillator at k > in H c s(t). The second 
is in the fixed basis in which (L|) = n^ — n _£- This 
difference is conserved as the scattering processes that 



lead to exchange of momenta between pairs of bosons is 
suppressed at infinite N. 

There are therefore half as many conserved quantities 
as the number of degrees of freedom. This is reminiscent 
of an integrable system in which the number of conserved 
quantities equals the number of degrees of freedom. For 
example, the free field theory (A = 0) is integrable and 
has two conserved quantities for each (k, —k) pair: the 
energy in the 2d oscillator and the angular momentum 
L|. In such systems, the correlations in the late time 
state are conjectured to be reproduced by the generalized 
Gibbs ensemble (GGE^H : 



late-time state has more memory than the GGE about 
the initial conditions and the dynamics does not relax to 
the GGE in Eq. (pi). 



This leaves open the possibility that there is a more 
restrictive GGE that includes some unknown additional 
conserved quantities that will describe the late time 
states. While we cannot rule this out, we (and previ- 
ous workers on this problem) have not been able to find 
such quantities. 



IV. QUENCHES TO THE ORDERED PHASE 



M 



Pgge = exp(- > fiiOi) 



where the Oi commute with the Hamiltonian and the p^ 
are fixed by (Oi) in the initial state. This conjecture has 
been checked for the two-point correlation functions in 
the free field theory (if Wick's theorem is satisfied in the 
initial state, higher point correlation functions would also 
be correctly reproduced). 

The collection of conserved quantities in the O(N) 
model at infinite N suggests that the late time state 
should thermalize to the following GGE: 



Pgge = exp 



-f3H 



fc>0 



hlI 



(19) 



We argue in two (possibly related) ways below that the 
dynamics fails to relax to this GGE. 

First, note that Eq. ( 16 1 depends on a c-number as- 



sociated with the initial state. A reasonable expectation 
would be that correlation functions at late times also de- 
pend on this c-number. However, correlation functions 
derived from the GGE do not. This suggests that the 
GGE is not the correct description of the late-time dy- 
namics. In Appendix [SI we sharpen this intuition by 
showing that the structure factor at late times depends 
on A/g. It follows that two-point function in real-space is 
also a function of the A/g, thus proving the claim that the 
system does not relax to the GGE for all local observ- 
ables. We note that Ref. [27] arrived at the same conclu- 
sion by constructing a particular non-linear combination 
of equal-time correlation functions that were independent 
of time and involved Mr. 

The second argument hinges on the numerical obser- 
vation that m 2 s (t) generically approaches a non- negative 
constant, ml, as t — > oo. This implies that the different 

k modes decouple and the theory is effectively free at 
late times. Thus, an extensive set of conserved quanti- 
ties (equal to the number in the GGE) emerges at late 
times: 



E% = 



nd 



(1*1 



»?)l*sl a 



(20) 



and the system appears integrable. These extra con- 
served quantities are independent of (L~) : they are func- 
tions of c-numbers like Mr in the initial state. Thus, the 



1. Problem, methods and background 

We now turn to the late time behavior of proto- 
cols starting from ground states in the disordered phase 
wherein the coupling is changed to a value in which the 
ground state is ordered (see Figs. [l]and[2|. For speci- 
ficity and because our interest here is primarily in the 
late time behavior, we will take these protocols to be 
sudden quenches, i.e. the coupling will be changed in- 
stantly to its new value. Our conclusions will generalize 
mutatis mutandis to more elaborate protocols such as 
those of interest in the Kibblc-Zurek problem. Right af- 
ter the quench, the system is in an excited state for the 
new Hamiltonian, i.e. one with a non-zero energy den- 
sity. For a system which thermalizes, the late time state 
can be read off from the equilibrium phase diagram by 
converting the energy density into an equivalent temper- 
ature. Interestingly, despite the lack of relaxation in the 
quantum O(N) model at infinite N, we will find that the 
phase diagram is still a good heuristic for the late time 
behavior when the system is quenched, in a sense we will 
shortly make precise. 

But first, let us describe the setup of the problem more 
precisely. Recall that the mode functions are constrained 
by Eqs. (12 1 and (13 1. Polar coordinates make these 
constraints more transparent. Define fr = \fr\exp(i6r). 



Eq. ( 12 ) then provides a relation between the amplitude 



and the phase of the mode function: 

n _ "I 
" \k?' 



(21) 



Together with the initial condition 6r(Q) = 0, this allows 
us to reconstruct the mode function at all t from |/e(t)| 
alone. The equation of motion for |/n(£)| follows from 
Eq. ( fl5fl6| >: 



dt 2 



+ \ k \'+mi s (t) |/g|- 



1 



I/a. 



m 2 cS {t) = r + 



X 



d d k 



(27T 



- d \h\ 2 ^- 



= (22) 



(23) 



It suffices to solve these equations subject to the initial 
conditions discussed below to determine fr(t). 

We pick the initial state of the system for t < to be 
the ground state in the disordered phase. That is, we 



prepare the system in the ground state of H (Eq. (|1| 
with bare coupling r = r$ greater than the critical cou- 
pling r c . This state is the ground state of a free field 
theory with the effective mass ttiq. The relation between 
mo and r$ is given by Eq. (HI). In this state: 



A/"r = i, n kQ = \\k\z 



hit < o) 



hit < o) = -i 



(24) 



f2fco above is the frequency of the harmonic oscillator at 
k. At t = 0, we suddenly quench to the ordered phase 
(r < r c ). We note that in terms of our discussion earlier, 
the conserved angular momenta are all zero so that the 
associated chemical potentials /ig in the GGE are all also 
zero. 

For the most part we will rely on numerical solutions 
of the dynamical equations. Specifically, we numerically 
solve Eqs. ( 22|23 ) subject to the initial conditions in 
Eq. (24). We sample |/^| on a grid of points in momen- 



tum space with the infra-red spacing 1/L and the ultra- 
violet cut-off A; we present data only for the largest sys- 
tem sizes L in which the finite-size effects are minimaP^. 
We have also ascertained that the average energy is con- 
served at least to one part in 10 4 . The dimensional de- 
pendence of the late time physics is easily accessible as 
the spatial dimension d can be varied continuously in 
these simulations. As we employ a lattice in momentum 
space, we drop the delta function factors in the contin- 
uum theory from this point on. Correlation functions 
that appear in the rest of the article should rightly be 
thought of as structure factors. 

As advertised, we are looking to use the equilibrium 
phase diagram to rationalize the behavior of our solu- 
tions. To this end, we quickly review the properties of 
the finite T equilibrium phases (see Figs. fT] and [2]). First, 
in d < 2, the system is always disordered for all T > 0. 
Correlations decay to zero exponentially with distance 
and the inverse of the effective mass, l/m c g plays the 
role of the correlation length. m c g > in this phase. In 
d > 2, there are two phases for r < r c (r c is the zero 
temperature quantum critical point). For T > T c (r), the 
system is disordered. As we just discussed, m e g > in 
this phase. For T < T c (r), the system is ordered. The 
magnetization (<£>) is non-zero and the two-point func- 
tion ($(x) • $(0)) decays as a power law l/|x| d ~ 2 to the 
positive constant K^)! 2 - A numerically useful property 
of such a two-point function is that the volume under 
the curve on long length scales |f| scales as ~ |(<&)| 2 |r| d . 
The effective mass squared is zero in this phase and is, 
physically, the mass of the (N — 1) Goldstone modes. At 
the critical temperature T c , m 2 ff = and the two-point 
function decays to zero as l/|a?| rf 2 . 



2. Results 

Our first result is the one that we alluded to in Sec. lIHI 
the late time effective mass squared tends to a non- 
negative constant, ml > 0. This is observed in all our 
numerical solutions. Certainly, ml < is ruled out on 
grounds of stability. Heuristically, the result follows from 
the averaging over many momentum modes which oscil- 
late at different frequencies. Once the effective mass set- 
tles, the density matrix in the eigenbasis of H dephascs 
and is effectively diagonal. 

The limiting behavior of the effective mass divides into 
two classes with very different physical content. 

In the first case, the effective mass squared tends to a 
positive constant {ml > 0). It is easily checked that this 
implies a finite correlation length at late times. The late 
time behavior is then qualitatively the same as that of the 
equilibrium disordered state. We emphasize that we are 
not implying thermalization: ml is not a function of the 
excess energy density and the angular momenta (LS) = 
alone. It depends on many properties of the initial state. 
Similarly the form of the late time correlations are not 
given by the appropriate GGE. 

In the second case, m 2 ff tends to zero as t — > oo and 
the correlation length diverges in the same limit. We 
shall refer to this behavior as coarsening. This case sub- 
sumes two cases — that of strict coarsening and that of 
coarsening to a critical state. By strict coarsening, we 
mean the analog in our system of the process in an er- 
godic system in which the system is quenched from the 
disordered to the ordered phase and the symmetry is only 
broken locally^. The domains of broken symmetry grow 
with time; their characteristic size at late times, l co (t), 
grows as t l ' Zd , where Zd is a dynamic exponent. Within 
each bubble, the two point function heals to its equilib- 
rium value |(<&)| 2 on a length scale £ -C l co {t)- On the 
longer length scale l co (t), it decays to zero. In contrast, 
by coarsening to a critical state, we mean the analog 
of the the system approaching a critical state with no 
length scale except l co (t) and hence no domains exhibit- 
ing equilibrium magnetization. In either case, the system 
is self-similar on the scale l co (t) and a dynamical scaling 
theory emerges in the limit t — > oo, \k\ — >• holding \k\ Zd t 
fixed. In this limit, the equal time structure factor has 
the scaling form: 



st, 



f 



<**(*)*£(*)> - W s g ^ 



l t), 



(25) 



where Q is the scaling function. When the system is 
strictly coarsening, the volume under the two-point func- 
tion should grow as l co {t) d ■ That is, the structure factor 
at zero momentum should as l co (t) d . Thus, S = d for 
strict coarsening. On the other hand, when the system is 
coarsening to a critical state, the growth of the volume 
under the two-point function with time is slower. Thus, 
5<d. 

To numerically confirm that the system is coarsening, 
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FIG. 3. Plot of rrv^git) vs t for a sudden quench from the 
disordered phase to the ordered phase in d = 1 + e. The green 
line shows m% GE , the effective mass squared predicted by the 
GGE. Inset shows late time. System parameters: e = 0.1, 
L = 1200, A = n, r = r c + 2, r = r c - 1, A = 1. 
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FIG. 4. Plots of ($jf(*)$t(£)) vs |jfc| at different times for 
a sudden quench from the disordered phase to the ordered 
phase in d — 1 + e. System parameters: e — 0.1, L = 1200, 
A = 7r, ro = r c + 2, r = r c — 1, A = 1. 



we will ask that the structure factor at late times have 
the scaling form predicted by Eq. (251. We will then 
use the value of 8 to differentiate between the two pos- 
sibilities. Let us now turn to the amplitude and spatial 
dimension dependence of the late time behavior of our 
sudden quenches. 

Consider first sudden quenches in d < 2. Relatively 
small system sizes are sufficient to determine the late- 
time behavior near d = 1 rather than at d = 2. We there- 
fore work at d = 1+e. Fig.kSphows m^ s (t) for d = 1.1; the 
inset contains the longest time behavior. Observe that 
m^ ff can be negative in the course of the evolution. This 
is not alarming because it is not the steady state behavior 
and does not imply instantaneous imaginary correlation 
lengths. The negative values of m^ s merely indicates 
that ($ 2 (t)) < (— r). The mode functions then grow 
exponentially, causing m^ ff to become positive. This ex- 
plains the initial oscillatory behavior of rn^ s . However, 
at late times, ml s settles to the positive value m 2 (see 
inset). The green line in both plots is the effective mass 



L GEE- 



The system def- 
is almost 



b GGE 



squared predicted by the GGE, 
initely does not relax to the GGE as 
two orders of magnitude larger than m ^ . The equal time 
structure factor is plotted in Fig.|4]at different times. By 
the Riemann-Lebesgue lemma, only the time-averaged 
structure factor contributes to the two-point function in 
real-space. Thus, the green curve is sufficient to under- 
stand the behavior of correlations in real-space. All three 
structure factors lead to exponential correlations in real- 



space. 

The positivity of m 2 f is not specific to our choice of 
initial conditions. Recall that stability arguments dictate 



that 



> if the effective mass goes to a constant 



7 
as t — > oo. 
physically impossible in d 



We now show analytically that to? 



< 2. 



is 
To this end, suppose 



/?? 



/ 



0. Then, the solution for ft(t): 



/, = A-cos(|fc|i) + 



r?- 

^sin(lfclt), 
1*1 



(26) 



where At and Br are complex valued functions that de- 
pend on the initial conditions. Whatever the detailed 
form, their amplitude must be finite and non-zero and 



k fc J 

to be consistent with 



Im[AtBt] = 1 (from Eq. (12)). The above solution has 



7 



determined through Eq. (16) 



(or Eq. (23)). For d < 2, this is impossible as the RHS 
is different from zero by an amount divergent in the long 
time limit: 



d d k 



(27T) 



'-iWXi-i 



2-d 



(27) 



Therefore, 



7 



> for d < 2 and the late time state is 



disordered. The argument at d = 2 is more delicate; it 
involves showing that the RHS is different from zero by 
a finite quantity. 



The above proof suggests that for d > 2 different 
quenches can lead to a vanishing mi and hence coars- 
ening. Indeed, our numerical results for d > 2 confirm 
this expectation and show two kinds of late time behav- 
iors. Shown in Fig. p5us the plot of m^ s (t) at late times 
for a sudden quench of "small amplitude" to the ordered 
phase in d = 3. We define a "small amplitude" quench to 
be a sudden quench in which the injected energy density 
is smaller than that in the critical ensemble at r. At late 
times, the effective mass is seen to oscillate about zero 
with a decaying amplitude. Thus, mi is indeed zero. 

In Fig.pj] we plot various structure factors: 1) at t = 0, 
2) the averaged structure factor at late times and 3) the 
connected structure factor predicted by the GGE. We 
compute the averaged structure factor only for the wave 
numbers |fc| > 0.1. For t > 25, their amplitudes, |/d, 
oscillate about the green curve. Again, by the Riemann- 
Lebesgue lemma, the green curve is sufficient to compute 
the late time two-point function in real-space at distances 
< 27r/0.1 in lattice units. Observe that despite the lack 
of relaxation to the GGE, the tail of the green curve falls 
off as l/|fc| 2 , just like the curve computed from the GGE 
(black). In the inset in Fig. 6l we plot |/o(t)| as a function 
of t. At early times, \fa(t) is order one as in the initial 
state. At late times, it grows linearly with t. This implies 
that the volume under the two point function |/ | 2 grows 
as t 2 . This is consistent with the scenario of coarsening 
to a non-equilibrium critical state, as we will show now. 

We now check if the structure factor has the scaling 
form predicted in Eq. ( 25 1 . Shown in Fig. [7] is the plot 
of the scaling function Q vs \k\t. The curves at different 

|fc| collapse when Zd = 1,6 = 2. As 6 < d, the system is 
coarsening to a non-equilibrium critical state, as opposed 
to strict coarsening (in agreement with Ref . f3T|) . Physi- 
cally, when Zd = 1, the size of a correlated region at time 
t is set by the horizon. This is the fastest possible growth 
and is possibly too fast to establish long-range order. An 
interesting avenue for future work is to investigate the 
nature of the coarsening process in the O(N) model at 
infinite N limit in different spatial dimensions. Prelim- 
inary results show that the system always coarsens to a 
non-equilibrium critical state. However, the critical state 
itself varies with dimension 37 . 

On injecting an energy density greater than that in 
the critical ensemble at r, the late time behavior is qual- 
itatively the same in Figs. |3|4| Again, as was the case 
for d < 2, the GGE does not reproduce the late time 
behavior. 



3. The step approximation 

The dimensional dependence of the late time physics is 
well captured by the step approximation, first introduced 
in Ref. [2H1 Here, one approximates m 2 ff (£) by a step 



0.015 



-0.015 




t 



FIG. 5. Plot of rrigg(i) vs t for a sudden quench of small 
amplitude (see text for definition) from the disordered phase 
to the ordered phase in d = 3. The ratio of the injected excess 
energy density to the critical energy density is approximately 
0.7. System parameters: L = 700, A — tt, tq — r c + 2, 
r = r c - 4, A = 1/2. 



function: 



'cir 



(*) 



'()• 



m~ 



t < 
t > 0. 



(28) 



The initial state fixes to 2 ,. This approximation is quite 
coarse in that it ignores the intricate early-time behavior 
of m 2 ff (t) visible in Fig.pl However, it builds in the late- 
time constancy of to 



cir * 



a key feature of the infinite 
N dynamics. Within this approximation, the equations 
of motion can be solved analytically as the dynamical 
problem is equivalent to a sudden quench in a free field 
theory, to 2 is a free parameter in this solution. It is fixed 



by requiring that the self-consistency relation in Eq. ( 23 ) 



hold at late times. For further details of the method, see 
Ref. [211 

Although m s is numerically not equal to the actual 
late time effective mass, to/, the two share qualitative 
features. Specifically, the dimensional dependence of toj 
discussed in the previous subsection is reflected by m s . 
First, we find that for d < 2, to 2 is always positive for 
any m . Thus, the final state is always disordered. Sec- 
ond, for d > 2, the value of to 2 depends on to 2 ,. For deep 
quenches or "large" too, m s > and the final state is 
disordered. On decreasing the amplitude of the quench, 
the value of m 2 decreases, until to 2 = 0. After this point, 
we find no solutions to the self-consistency relation and 
the approximation breaks down. This breakdown is in- 
dicative of the new physics of coarsening at late times. 
Note though that coarsening cannot be captured within 
this approximation as the precise way in which m^ s (t) 
approaches zero is important. 

To summarize, the step approximation leads to a dis- 
ordered state in d < 2 for any quench and in d > 2 for a 
deep quench. For shallow quenches in d > 2, the approx- 
imation breaks down, indicating new late-time behavior. 
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FIG. 6. Plots of ($g(t)$t(t)) vs |fc| at different times for a 
quench of small amplitude from the disordered phase to the 
ordered phase in d = 3. The green curve is the averaged 
structure factor for t > 25 (see text). Inset: The amplitude 
of the mode function at |fc| = 0, |/o(t)|, vs t. See the label of 
Fig. [5]for system parameters. 




FIG. 7. The structure factor multiplied by \k\ 2 vs \k\t for 
the quench of small amplitude from the disordered phase to 
the ordered phase in d — 3. The six curves that exhibit the 



scaling collapse predicted by Eq. ( 25 1 are at the six smallest 
values of |fc|, |fc| = 2-Km/L, where m — 1, ... 6 and L = 700. 
We conclude that za — 1> 8 — 2. 



SCALING AND QUENCHES TO THE 
CRITICAL POINT 



Thus far we have focused on the dimensional depen- 
dence of the late time physics upon quenching to the or- 
dered phase. Now we turn to the question of scaling for 
such quenches due to proximity to the critical point sepa- 
rating the ordered and disordered phases. This aspect of 
the physics has an interesting dimensional dependence of 
its own which is the non-equilibrium analog of the varia- 
tion of equilibrium critical behavior with dimension. We 
will specifically be interested in the reversion to Gaussian 
critical behavior in d > 3. We will focus almost entirely 
on quenches to the critical point and note at the end the 
generalization to quenches into the ordered phase. 

As before, we prepare the system in the ground state 
in the disordered phase at some r$ > r c . The initial 
conditions are once again as detailed in Eq. (24) 



At 
t = 0, we quench the system to the critical point. 

A first observation is that in any dimension the late 
time state is disordered, i.e. m? > for an initial non- 
zero rriQ. This is observed in our numerical solutions 
but it can also be argued, on analytic grounds parallel- 
ing our demonstration that coarsening is impossible upon 
quenching into the ordered phase. For d < 2, our pre- 
vious argument already suffices and for d > 2, one can 
extend it via a consideration of a larger set of momenta^ 
at the critical coupling. We note that this result is again 
consistent with reasoning based on the equilibrium phase 
diagram — injecting a non-zero energy density at the crit- 
ical coupling should lead to a finite correlation length. 
However, this result hides an important difference be- 
tween d < 3 and d > 3 to which we now turn. This 
difference is most sharply visible if we consider a scaling 
limit for sudden quenches, formulated in Ref. 14 which 
makes their universal physics manifest. 

Let us first consider 1 < d < 3 where the critical point 
is interacting even in the infinite N theory. When the 
initial state is in the vicinity of this critical point, its cor- 
relation length, l/m ' , diverges as (r — r c )~ v , where the 
critical exponents are v = l/(d — l),z = 1 at infinite N. 
On the length/time scale l/m 0l the equilibrium physics 
is universal. On the same scale, the dynamics following a 
sudden quench is also universal in an appropriate scaling 
limit. This is the limit of tq —$■ r c when lengths/times are 
measured in units of the initial correlation length/time, 
1/mo. In this limit, the one and two point functions of 
an operator O have the scaling forms: 



(O(x,t)) mo 
(O(f,t)0(0,0)) mo 



m g o (tm ) 
ml A g oo (\x\m a ,tm ), 



(29) 
(30) 



where A is the scaling dimension of the operator O and 
we have used the translational and rotational invariancc 
of the system. The out-of-equilibrium physics in this 
limit is conjectured to be universal because the quench 
is very shallow - only the universal low-energy long- 
wavelength part of the energy spectrum at the critical 
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FIG. 8. Plots of mlff/iriQ vs tmo a t different values of L% 
in d=2. At each Lmo, we show the behavior for three (small) 
values o f m p. The scaling collapse provides strong evidence 
for Eq. (311. System parameters: A 



1, r 
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FIG. 9. M{oo) vs spatial dimension d computed for m^ 1 = 
20, 40, 80 with Lm — 10. T he soli d line is M(oo) obtained in 
the step approximation (Sec. |IV"3| ) at the same system system 
sizes and initial conditions. System parameters: A = 1, r = 
r c , A is d-dependent and chosen to ensure that the initial state 
is in the equilibrium scaling region. 



point is excited for t > 0. Those readers familiar with 
the Kibble-Zurek mechanism will notice that the above 
scaling limit is very similar to the one in slow ramps. The 
correlation length and time in the initial state play the 
roles of the Kibblc-Zurck length and time. 

The previous computations of universal correlation 
functions in d = 1 are in the scaling limit described 
above. For example, the two-point corre lation functions 
computed by Calabrese and Cardjff^H us i n g boundary 
conformal field theory have the scaling form in Eq. ( 30 1 . 



Finally, the formalism above can be easily generalized to 
the cases when (a) the dynamic exponent z ^ 1, (b) when 
the system size is finite and (c) when the quench is not 
exactly to the critical point. 

To test the scaling hypothesis, we consider the one- 
point function, (\<f>\ 2 (t)) mo / N . Equivalently, consider the 



effective mass squared m^ ff . Accounting for the finite 
system size L in the numerics, the scaling form of m^g 
is: 



ml s (t,L;mo) ~ mlM(tm ,Lm ). 



(31) 



We discuss the early and late time behavior of M. below. 
By construction in the thermodynamic limit, A4 (0~) = 
1. At t = + , m^ s is negative. This reflects the finite 
correlations in the initial state, unlike the ground state at 
the critical point. Using the form of /e(0~) in Eq. (23 1, 
it is easy to see that: 

M(Q + ) ~ -m*- 3 . (32) 

Thus, .M(0 + ) is negative and divergent in the scaling 
limit. 

The full scaling function can be determined numeri- 
cally. In Fig. [8j we show M. for Lmo = 6, 20 in d = 2. 
The data collapses for three different values of mo at each 
Lmo, providing strong evidence for Eq. (31 1. The nega- 



tive divergence of M. near + is not shown in order that 
the late time behavior be visible. Generically, we find 
that M. tends to a positive constant at late times. Thus, 
the late time state is indeed disordered. 

In Fig. |j] M(oo) is plotted as a function of dimen- 
sion for different values of mo at fixed Lmo — 10. For 
all d < 3, the extrapolation to the thermodynamic limit 
leads to a positive value of A4 as tm —> 00. The solid 
line is .M(oo) in the step approximation discussed in 



Sec. IV 3 In agreement with Eq. (31), .M(oo) is finite 



in the dynamic scaling limit. The close agreement with 
the numerical data points appears fortunate rather than 
principled. 

Now consider the case when d > 3 where the fixed 
point is Gaussian. If we work directly at the fixed point, 
we have a free field theory, and the time-dependent mode 
functions can be calculated exactly for the sudden quench 
to the critical point (see, for example Ref. 132]) . In the 
language of this article, the computation amounts to ig- 
noring the self-consistency equation for the effective mass 
and assuming that it is zero at all times. The result is 
that the correlation functions at late times have a power- 
law behavior in real space at long distances (light cone 
effects are irrelevant to this discussion). Now, this is at 
odds with the observation of a finite correlation length 
in Fig. [9] at any non-zero mo- This discrepancy is re- 
solved if we take proper account of the irrelevant quar- 
tic term which is dangerously irrelevant already for the 
equilibrium behavior, i.e. it cannot be neglected to get 
proper asymptotic results. The same is true for the non- 
equilibrium dynamics. This is easy to see in the step 
approximation. The dependence of the late time (finite) 
correlation length, £/, on the initial correlation length 
l/m is given by: 



£/«i 



-(d-3)/2 



as mo — > 0. 



(33) 



In the marginal d — 3 case, £/ is only logarithmically 
greater than mo- Thus, if we take the scaling limit keep- 
ing distances on the order of 1/mo fixed as described in 
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Ref. Q31 we will indeed find Gaussian behavior but this 
will no longer be a good description of the truly long time 
asymptotics on the longest time and length scales. As in 
typical problems with dangerously irrelevant variables, 
there are now two divergent scales to contend with — a 
phenomenon likely to be much more common in the non- 
equilibrium setting as has been noted already with a dif- 
ferent example involving string-net coarsening in Ref. 1381 
Finally, we note that the dangerous irrelevance of in- 
teractions is also germane to quenches into the ordered 
phase where again a scaling limit can be defined as above. 
Again, we will find that in d > 3 Gaussian results hold 
in the scaling regime, which will exhibit an exponential 
growth of the local order parameter before being cutoff at 
the parametrically longer scale by the coarsening physics 
we described previously. 



VI. CONCLUDING REMARKS 

In our current understanding, as sketched in this pa- 
per, the quantum 0(N — oo) vector model appears in- 
termediate between generic systems that exhibit thermal- 
ization starting out of equilibrium and integrable systems 
that do not. It does exhibit stationary behavior at long 
times following parameter changes but it does not exhibit 
thermalization. This behavior is consistent with it ap- 
pearing to have only half the number of conserved quan- 
tities appropriate for a fully integrable system. Interest- 
ingly, the GGE constructed from the known conserved 
quantities does not describe the late time stationary 
states which is inconsistent with the GGE conjectures'. 
A definitive resolution to the question of whether the 
GGE conjecture is false or whether there are additional 
conserved quantities is desirable. In the absence of the 
latter the 0(N = oo) vector model would provide an 
example of a system exhibiting a large number of emer- 
gent conserved quantities at late times. In such a case, 
at least the most salient properties of the system might 
be captured in terms of a modified GGE which includes 
chemical potentials for the emergent conserved quanti- 
ties. 

Interestingly, the late time states of the 0(N = oo) vec- 
tor model following quenches exhibit "knowledge" of its 
equilibrium phase diagram. They exhibit coarsening if 
and only if the model exhibits a finite temperature or- 
dered phase, which in turn depends on the dimensionality 
of the system. However, this coarsening process appears 
to be to a state that is critical. That is, it does not show 
any signatures of developing long range order on length 
scales smaller than the coarsening length scale. Elucidat- 
ing the nature of this critical coarsening and its precise 
dependence on initial conditions is a fit subject for future 
work. 

We note that the temporal structure of relaxation ap- 
pears to be quite complicated at large but finite N. The 
infinite N theory already exhibits finite times scales for 
the appearance of stationary states in the absence of 



coarsening and algebraic relaxation in the presence of 
coarsening. This behavior is initial state dependent. For 
large N we expect this behavior to give way to genuine 
thermalization, at least for d > 1, on a time scale that 
is parametrically large in N. How much of this intricate 
time dependence survives to, say, N = 3 is an interesting 
question. 

Finally we note that in d — 3 the critical fixed point 
is Gaussian. In its neighborhood, the scaling limit will 
be described by the Gaussian theory as we illustrated in 
this paper for quenches to the critical point. This is po- 
tentially amenable to experimental work in cold atomic 
systems where the N = 2 case is the critical theory of 
the Mott insulator to superfluid transition with particle 
hole symmetry. 
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Appendix A: Lack of relaxation to the GGE 

In this appendix, we show that the correlations in the 
late time state of the O(N) model are not reproduced 
by the appropriate GGE Eq. ( |19[ ). Define a new set 
of variables linearly related to the mode functions ft in 
Eq. (fill: 



(Al) 

We interpret g^ as the complex coordinate of a 2d classi- 
cal particle labelled by k > 0. Next, define the momen- 
tum of each particle to be: 



■k = g 



(A2) 



When absent, the momentum subscript is implied. Then, 
the condition in Eq. (12 1 implies that the angular mo- 



mentum of the particle at k > is given by: 
I* = -Im[g(t)n(t)} = -AT. 



(A3) 



This classical angular momentum should not be confused 
with the conserved angular momentum L z in Eq. (17 1. 



They are completely unrelated. The above equation im- 
plies that l z is a constant of the motion for the classical 
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system and is set by the value of Af. The dynamics of 
the classical particles is governed by the Hamiltonian: 



H r< 



\A 2 + {\k\ 2 +r)\g\ 2 + 



A 



fc>0 



q>0 



(A4) 



It is easily checked that the equations of motion derived 
from H c i reproduce Eq. ( 15|16 ). On solving for the dy- 
namics of the classical particles generated by H c i with 
fixed angular momentum V and the initial conditions 
derived from Eq. ( |14[ ), we can compute any observable 
in the O(N) model through Eq. (11). This classical in- 
terpretation is useful because it unambiguously identifies 
the role of the c-number Af in the solution. The radial 



coordinate of the classical particle, |g(£)|, depends on the 
conserved angular momentum l z at late times. As the 
equal time structure factor S(k, t) is given by: 



S(M) 



where ($ r $L) 
\ k k>i 



lfls(*)l ! 



S(k,t)(2ir) d S d (k-k') 



(A5) 
(A6) 



the time-averaged structure factor at late times must de- 
pend on ZS as well. The structure factor in the GGE 
however is independent of this quantity. Thus, the two 
do not agree and the dynamics in the O(N) model does 
not thermalize in the generalized sense. 
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